Single-cell profiling of transcriptome and histone modifications with EpiDamID.

Recent advances in single-cell sequencing technologies have enabled simultaneous measurement of multiple cellular modalities, but the combined detection of histone post-translational modifications and transcription at single-cell resolution has remained limited. Here, we introduce EpiDamID, an experimental approach to target a diverse set of chromatin types by leveraging the binding specificities of single-chain variable fragment antibodies, engineered chromatin reader domains, and endogenous chromatin-binding proteins. Using these, we render the DamID technology compatible with the genome-wide identification of histone post-translational modifications. Importantly, this includes the possibility to jointly measure chromatin marks and transcription at the single-cell level. We use EpiDamID to profile single-cell Polycomb occupancy in mouse embryoid bodies and provide evidence for hierarchical gene regulatory networks. In addition, we map H3K9me3 in early zebrafish embryogenesis, and detect striking heterochromatic regions specific to notochord. Overall, EpiDamID is a new addition to a vast toolbox to study chromatin states during dynamic cellular processes.

Simultaneous quantification of protein-DNA interactions and transcriptomes in single cells with scDam&T-seq.

Protein-DNA interactions are essential for establishing cell type-specific chromatin architecture and gene expression. We recently developed scDam&T-seq, a multi-omics method that can simultaneously quantify protein-DNA interactions and the transcriptome in single cells. The method effectively combines two existing methods: DNA adenine methyltransferase identification (DamID) and CEL-Seq2. DamID works through the tethering of a protein of interest (POI) to the Escherichia coli DNA adenine methyltransferase (Dam). Upon expression of this fusion protein, DNA in proximity to the POI is methylated by Dam and can be selectively digested and amplified. CEL-Seq2, in contrast, makes use of poly-dT primers to reverse transcribe mRNA, followed by linear amplification through in vitro transcription. scDam&T-seq is the first technique capable of providing a combined readout of protein-DNA contact and transcription from single-cell samples. Once suitable cell lines have been established, the protocol can be completed in 5 d, with a throughput of hundreds to thousands of cells. The processing of raw sequencing data takes an additional 1-2 d. Our method can be used to understand the transcriptional changes a cell undergoes upon the DNA binding of a POI. It can be performed in any laboratory with access to FACS, robotic and high-throughput-sequencing facilities.

Simultaneous quantification of protein-DNA contacts and transcriptomes in single cells.

Protein-DNA interactions are critical to the regulation of gene expression, but it remains challenging to define how cell-to-cell heterogeneity in protein-DNA binding influences gene expression variability. Here we report a method for the simultaneous quantification of protein-DNA contacts by combining single-cell DNA adenine methyltransferase identification (DamID) with messenger RNA sequencing of the same cell (scDam&T-seq). We apply scDam&T-seq to reveal how genome-lamina contacts or chromatin accessibility correlate with gene expression in individual cells. Furthermore, we provide single-cell genome-wide interaction data on a polycomb-group protein, RING1B, and the associated transcriptome. Our results show that scDam&T-seq is sensitive enough to distinguish mouse embryonic stem cells cultured under different conditions and their different chromatin landscapes. Our method will enable the analysis of protein-mediated mechanisms that regulate cell-type-specific transcriptional programs in heterogeneous tissues.

Genome-lamina interactions are established de novo in the early mouse embryo.

In mammals, the emergence of totipotency after fertilization involves extensive rearrangements of the spatial positioning of the genome1,2. However, the contribution of spatial genome organization to the regulation of developmental programs is unclear3. Here we generate high-resolution maps of genomic interactions with the nuclear lamina (a filamentous meshwork that lines the inner nuclear membrane) in mouse pre-implantation embryos. We reveal that nuclear organization is not inherited from the maternal germline but is instead established de novo shortly after fertilization. The two parental genomes establish lamina-associated domains (LADs)4 with different features that converge after the 8-cell stage. We find that the mechanism of LAD establishment is unrelated to DNA replication. Instead, we show that paternal LAD formation in zygotes is prevented by ectopic expression of Kdm5b, which suggests that LAD establishment may be dependent on remodelling of H3K4 methylation. Our data suggest a step-wise assembly model whereby early LAD formation precedes consolidation of topologically associating domains.

Massive reshaping of genome-nuclear lamina interactions during oncogene-induced senescence.

Cellular senescence is a mechanism that virtually irreversibly suppresses the proliferative capacity of cells in response to various stress signals. This includes the expression of activated oncogenes, which causes Oncogene-Induced Senescence (OIS). A body of evidence points to the involvement in OIS of chromatin reorganization, including the formation of senescence-associated heterochromatic foci (SAHF). The nuclear lamina (NL) is an important contributor to genome organization and has been implicated in cellular senescence and organismal aging. It interacts with multiple regions of the genome called lamina-associated domains (LADs). Some LADs are cell-type specific, whereas others are conserved between cell types and are referred to as constitutive LADs (cLADs). Here, we used DamID to investigate the changes in genome-NL interactions in a model of OIS triggered by the expression of the common BRAFV600E oncogene. We found that OIS cells lose most of their cLADS, suggesting the loss of a specific mechanism that targets cLADs to the NL. In addition, multiple genes relocated to the NL. Unexpectedly, they were not repressed, implying the abrogation of the repressive activity of the NL during OIS. Finally, OIS cells displayed an increased association of telomeres with the NL. Our study reveals that senescent cells acquire a new type of LAD organization and suggests the existence of as yet unknown mechanisms that tether cLADs to the NL and repress gene expression at the NL.

Genome-wide maps of nuclear lamina interactions in single human cells.

Mammalian interphase chromosomes interact with the nuclear lamina (NL) through hundreds of large lamina-associated domains (LADs). We report a method to map NL contacts genome-wide in single human cells. Analysis of nearly 400 maps reveals a core architecture consisting of gene-poor LADs that contact the NL with high cell-to-cell consistency, interspersed by LADs with more variable NL interactions. The variable contacts tend to be cell-type specific and are more sensitive to changes in genome ploidy than the consistent contacts. Single-cell maps indicate that NL contacts involve multivalent interactions over hundreds of kilobases. Moreover, we observe extensive intra-chromosomal coordination of NL contacts, even over tens of megabases. Such coordinated loci exhibit preferential interactions as detected by Hi-C. Finally, the consistency of NL contacts is inversely linked to gene activity in single cells and correlates positively with the heterochromatic histone modification H3K9me3. These results highlight fundamental principles of single-cell chromatin organization. VIDEO ABSTRACT.

Single-cell dynamics of genome-nuclear lamina interactions.

The nuclear lamina (NL) interacts with hundreds of large genomic regions termed lamina associated domains (LADs). The dynamics of these interactions and the relation to epigenetic modifications are poorly understood. We visualized the fate of LADs in single cells using a "molecular contact memory" approach. In each nucleus, only ~30% of LADs are positioned at the periphery; these LADs are in intermittent molecular contact with the NL but remain constrained to the periphery. Upon mitosis, LAD positioning is not detectably inherited but instead is stochastically reshuffled. Contact of individual LADs with the NL is linked to transcriptional repression and H3K9 dimethylation in single cells. Furthermore, we identify the H3K9 methyltransferase G9a as a regulator of NL contacts. Collectively, these results highlight principles of the dynamic spatial architecture of chromosomes in relation to gene regulation.

Role for Msh5 in the regulation of Ig class switch recombination.

Ig class switch recombination (CSR) and somatic hypermutation serve to diversify antibody responses and are orchestrated by the activity of activation-induced cytidine deaminase and many proteins involved in DNA repair and genome surveillance. Msh5, a gene encoded in the central MHC class III region, and its obligate heterodimerization partner Msh4 have a critical role in regulating meiotic homologous recombination and have not been implicated in CSR. Here, we show that MRL/lpr mice carrying a congenic H-2(b/b) MHC interval exhibit several abnormalities regarding CSR, including a profound deficiency of IgG3 in most mice and long microhomologies at Ig switch (S) joints. We found that Msh5 is expressed at low levels on the H-2(b) haplotype and, importantly, a similar long S joint microhomology phenotype was observed in both Msh5 and Msh4-null mice. We also present evidence that genetic variation in MSH5 is associated with IgA deficiency and common variable immune deficiency (CVID) in humans. One of the human MSH5 alleles identified contains two nonsynonymous polymorphisms, and the variant protein encoded by this allele shows impaired binding to MSH4. Similar to the mice, Ig S joints from CVID and IgA deficiency patients carrying disease-associated MSH5 alleles show increased donor/acceptor microhomology, involving pentameric DNA repeat sequences and lower mutation rates than controls. Our findings suggest that Msh4/5 heterodimers contribute to CSR and support a model whereby Msh4/5 promotes the resolution of DNA breaks with low or no terminal microhomology by a classical nonhomologous end-joining mechanism while possibly suppressing an alternative microhomology-mediated pathway.